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Abstract 

We prepese seme methods for studying hadronic sequential two-body decays 
involving more spinning particles. It relies on the analysis of T-odd and T-even 
asymmetries, which are related to interference terms. The latter asymmetries turn 

out to be as useful as the former ones in inferring time reversal violating 
observables; these in turn may be sensitive, under some particular conditions, to 
possible contributions beyond the standard model. Our main result is that one can 
extract such observables even after integrating the differential decay width over 
almost all of the available angles. Moreover we find that the correlations based 
exclusively on momenta are quite general, since they provide as much information 
as those involving one or more spins. We generalize some methods already proposed 
in the literature for particular decay channels, but we also pick out a new kind of 
time reversal violating observables. Our analysis could be applied, for example, to 

data of LHCb experiment. 
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1 Introduction 



Interference is a typical quantum mechanical effect and can be exploited experimen- 
tally to detect important phenomena. For example, parity violation was discovered 
in the fifties by hypothesizing[l] , and successively exhibiting [2], an interference term 
between two decay amplitudes which behave differently under parity inversion; these 
were identified later [3] as matrix elements of the vector and axial current respectively. 
As is well-known, such a term is proportional to a scalar product of the type J ■ p[l], 
where J and p are respectively an angular momentum and a momentum of parti- 
cles involved in the process under study. It is worth recalling that this dependence 
was picked out by a suitable asymmetry [2]. Similarly, the asymmetries connected to 
direct CP violation — revealed, it is worth mentioning, only many years after the 
discovery of indirect CP violation in the Kl — > tttt decays [4, 5] — can be detected 
if two different amplitudes, with different weak and strong phases, contribute to the 
decay [6, 7]. 

It could be tempting to propose for time reversal violation (TRV) something simi- 
lar to parity violation. That is, wc could study interference terms between amplitudes 
which behave differently under time reversal (TR). These terms correspond to corre- 
lations of the type[8] 

C = Vi • V2 X V3, (1) 

where Vj [i = 1,2,3) are either momenta or angular momenta of particles involved 
in the process. Such correlations, named T-odd[7, 9, 10, 11, 12], can be defined in 
processes where more than two particles are involved. Correlation (1) may be revealed 
by the asymmetry [8] 

^ {C>0)-{C< 0) 

(C>0) + (C<0)" ^ ^ 

However in this case we are not so lucky as with parity violation. Indeed, the genuine 

TRV effects are mimicked by fictitious T-odd interference terms, caused for example 

by strong or electromagnetic spin-orbit interactions[13, 14]. Such effects, present in 
scattering proccsscsflO, 11, 12, 13, 14] as well as in decays (see for example refs. 
[7. 8, 9, 15, 16, 17, 18, 19, 20, 21, 22, 23]), produce the same T-odd correlations 
(1) as real TRV. In particular, they mask real TRV in a weak decay, owing to final 
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state interactions (FSI). It is therefore quite a serious problem to separate the two 
contributions [7, 8, 24]. But combining the asymmetry of a given triple product with 
the one of its CP-conjugated gives rise to a real TRV[8, 24, 25]. 

CP violations (CPV), as well as TRV, which can be seen as the counterpart of 
CPV owing to the CPT symmetry, are often used for investigating possible clues 
of new physics (NP). Indeed, baryon-antibaryon asymmetry in the Universe cannot 
be explained by the sole CKM mechanism of CPV[26]. A common strategy to find 
hints beyond the standard model (SM) consists of investigating experimentally those 
processes where the CPV (or TRV) are tiny in the SM predictions. Indeed some of 
such violations may indicate deviations from the SM. Especially, we recall B — > nK 
decays[27, 28, 29], Bs — Bs mixing[30, 31] and the like-sign dimuon asymmetry [32]. 
However till now experimental and theoretical uncertainties do not allow any conclu- 
sive results. Confirmations to the first few discrepancies are demanded. In this sense, 
experiments on sequential two-body decays, mainly of the type 

^ V^V2 (3) 

(with Vi^2 vector mesons), have been performed[17, 28, 30, 31] and suggested[21, 33, 
34, 35, 36, 37, 38, 39]. Such types of decays, where more spinning particles are 
involved, offer the advantage of the angular correlations [7, 8, 24, 25, 40, 41, 42, 43], 
which may provide the moduli and the relative phases of the helicity or transversity 
amplitudes [44, 45]. It is worth noting, in the context of such decays, that the above 
mentioned T-odd correlations are more and more frequently proposed[18, 19, 20, 21, 
22] and also used[15, 16, 17] for discovering possible clues of NP. This is because 
these correlations may be preferred to direct CP-asymmetries in the cases when the 
relative strong phases are neghgibly small[25, 46]. Moreover, as regards the helicity 
amplitudes, they are a further, suitable tool for testing the SM[38, 39, 47], possibly 
by examining fake [38, 48] T-odd observables. 

The aim of the present paper is to illustrate some methods for analyzing two- 
body hadronic weak decays. In particular we study the properties of asymmetries 
connected to T-odd products relative to a sequential two-body decay; moreover, we 
consider some T-even asymmetries, seldom used in the literature but equally useful. 
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First of all, we define suitable T-odd and T-even correlations, of the type (1) or 
slightly more complicated, and the corresponding asymmetries. Then we show how 
to use such asymmetries for extracting observables sensitive to TRV and, possibly, 
to NP. In particular, we illustrate some tests for probing clues to physics beyond the 
SM. This work implements in some way a previous one[49], where we proposed tests 
based on the knowledge of moduli and relative phases of the amplitudes of decays of 
the type considered here. Moreover it generalizes some of the methods adopted in 
the hterature in decays of the type (3). 

The paper is organized as follows. Sects. 2 to 4 are devoted to the definitions and 
to the analytical expressions of single and double T-odd and T-even asymmetries, 
written as functions of helicity amplitudes. Several TRV observables are defined in 
sect. 5. In sect. 6 we carry out an explicit example. Sect. 7 is dedicated to relations 
with previous papers. Various tests of the SM are presented in sect. 8. Lastly some 
conclusions are drawn in sect. 9. 



2 Single T-odd and T-even Asymmetries 

Here and in the following section we illustrate some T-odd and T-even asymmetries, 
which can be inferred from a sequential decay of the type 

J — >■ a 6, a — >■ ai a2, 6 — >■ 6i 62- (4) 

Such observables depend on correlations of the type (1) or similar. The T-odd asym- 
metries have been proposed by several authors, in different contexts. The most re- 
cent contributions regard hadronic[21, 23, 24, 25, 34, 38, 40, 41, 50, 51, 52] and 
semi-leptonic[53, 54, 55, 56, 57] decays, searches for top decays[9, 58, 59] and for 
new particles[18, 20, 22, 60]. Also some experiments[15, 16, 17] adopt this tech- 
nique. On the contrary, T-cvcn asymmetries are not so frequently suggested in the 
literaturc[9, 51, 61]. In order to introduce these observables, we define preliminarily 
suitable reference frames. 
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2.1 Reference Frames 



First of all, we define a canonical frame, at rest with respect to the parent resonance 
J: 

Pin ^ Pin ^ PJ ^ ^ ^ /_s 

y=i r, z=i r, x = yxz, (5 

|Pm| I Pin X PjI 

where pj„ and pj arc, respectively, the momenta of the initial beam and of the 
resonance J in the laboratory frame. 

On the other hand, the successive decays of the particles a and b are more conve- 
niently described in the helicity frames. For example, as regards particle a, it is given 
by the following three mutually orthogonal unit vectors: 



Gl = 1 — r, Gt = -p —r, ejv = er X e^, (6) 



Pa 






er = 77 7—r 


IPal' 


z X 



Pa being the momentum of a in the canonical frame. Analogous definitions hold true 
for particle b. 



2.2 Definition of Single Asymmetries 
2.2.1 Single T-odd Asymmetries 

Consider the scalar product 

7f = Pa, • e^, (7) 

where p^j is the momentum of the particle Oi in the rest frame of a. This is a T-odd 
quantity. Correspondingly, we define the T-odd asymmetry 

^ ^ iV(Ti^ > 0) - iV(Ti^ < 0) 
" N{T^ > 0) + N{T^ < 0) ' ^ ' 

where N{T^ > 0) [N{T^ < 0)] is the number of decays such that, for a given p^^, 
the scalar product above is positive (negative). Another, independent, T-odd product 
can be defined as 

T,^ = p,,.e;v, (9) 

where Pf,^ is the momentum of particle bi in the rest frame of b. Obviously, this 
product induces an asymmetry analogous to (8). 
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If different than zero, such asymmetries do not imply TRV, since, as told, they may 
derive contributions also from strong or/and electromagnetic interactions [7, 8, 13, 14]. 
However, as we shall see in sect. 5, it is possible to define quantities sensitive to such 
a violation[8, 24, 25]. 

2.2.2 Single T-even Asymmetries 

The scalar products 

rj(^) = Pa, • BTiD (10) 

are T-even quantities. Analogously to the T-odd asymmetries of subsect. 2.2.1, we 
define the T-even asymmetries 

^T(L) ^ NiT^L) > 0) _ Ar(Tj(^) < 0) 
7V(7f ^''^ > 0) + 7V(7f < 0) ■ 

Also in this case it is possible to define two more (T-even) asymmetries by substituting 
Pb^ to Pa^ in the definitions (10). 



2.3 Differential Two-Body Decay Width 

We deduce, here and in the following subsections, the expressions of the differential 
width of a two-body decay and of the asymmetries just defined before. The calcula- 
tions are based on the formalism of the density matrix[62, 63, 64, 65]. The starting 
point is the differential decay width for sequential two-body decays [66, 67]. It reads 
as 



1 d^T 



aia2bib2 



= NjabWY.<<%a,.{^^^a.^b). (12) 



IJ-a,IJ-b 



Here the F's are the partial decay widths of the decays (4); moreover 

2s. + 12s5 + 12J+l 

^ 4^^' ^^^^ 

W = \B{pmBipl)\'\Bipl)\', (14) 
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M,M' 

<A. = <A.(E KA.n-^/^ (16) 

< = E l^u" u l'( E l^u" ^ H"', (17) 

^^6 = l^^bjM6j^( 51 l^^biMfcJ^) (18) 

J, Sa and Sb are the spins, respectively, of the parent resonance and of particles a and 
b. ^AaAfc) ■^th.iHa2 ^^(-i^iba rotationally invariant two-body decay amphtudes 

relative to the three decays (4). They depend, respectively, on the hehcities Aa(A^) 
and Ab(A^), /lai and fj,a2, and ii^^ and /ib^, such that 

Aa - Ab = A, - A^ = A', - IJ'a2 = /^a, /^6i " I^b2 = /^6- (19) 

M(M') is the third component of the spin of J in the canonical frame. Q, Qa, 
denote the directions of, respectively, particles a, ai and hi in the rest frames of their 
parent resonances: as usual, we have set Q = {6, 0), where and (j) are, respectively, 
the polar and azimuthal angle; similar definitions hold for and 0.^. p^mm' 
B(jpys are, respectively, the spin density matrix of the parent resonance and the 
relativistic Breit-Wigner functions of the resonances, normalized as 

E&iP'y)-^, rdpW)?-i- (20) 

Lastly, T>'lj^{^l) is a Wigner "D- matrix element, defined as[68] 

-D'mA^) = e-'^^^di^iO), (21) 
di^{9) = {JM\e-^'y'\JA), (22) 

where Jy is the y-component of the angular momentum. 

Problems arise if one takes into account the energy dependence of the decay am- 
plitudes and of the density matrix; however, if the resonances are sufficiently narrow, 
as it usually happens, such a dependence can be neglected. Therefore, in an analysis. 
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it would be convenient to integrate over the energies of the resonances. But we do not 
perform such integrations, in order to point out an important theoretical subtlety, to 
be clarified in sect. 5. 

Now we integrate eq. (12) over d^Qad^^b- To this end, we have to take into 
account the normalization 

J d'nvi:^{n)vi,^{n) = ^^Smm' (23) 

for Wigner's "D-matrices. As a result we get the expression of the differential two-body 
decay width, i. e., 



Xa,\ M,M' 

Here 

.r 2J + 1 

Afj = (25) 

Rigorously speaking, the distribution F depends also on pj, pi and pi, through the 
factor W. However, in order to simplify the notation, we have omitted such a depen- 
dence in the argument of the distribution. Wc shall adopt this convention for all of 
the asymmetries that we shall define in the present paper. 

More important, we remark that the differential decay width (12), and therefore all 
of the expressions derived from it, do not depend strictly on the decay amplitudes, but 
rather on the dimensionless quantities (16); as we shall see, these play an important 
role in the definitions of the asymmetries and, from now on, will be named "reduced" 
amplitudes. 

2.4 Analytical Expression of the T-odd Asymmetry 

The expression of the observable r(Q)A^ (Q) reads 

Jo '-J-Tr/2 J-Kjl J J 

= •AA.a^E^EEeilA^.'^"^"" 



X <xyxtx, E K\{n)vU,{n)p%,. (26) 



M,M' 
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Here Ar(fi, fla^b) is a short notation for the differential decay width (12). Moreover 
we have set 

A/-.. = D. = {X^-K-l)/2, (27) 

e^x;,. = /J <».ii'«».<:.. (28) 

the d-matrices being defined by eq. (22). Lastly, we have denoted by Yl' the sum over 
those for which Xa — X'^ is odd. 

2.5 Expressions of the T-even Asymmetries 

The expression of the asymmetry is given by 



r(Q)A^(Q) = j\ineadea[f^ - £^]dct>aj d''n^/^v{n,nA) 



X 



Ma A„,A6 A^ ^« 

<xA, E ^^ma(^)^^m'A'(^)Pmm'- 

M,M' 



(29) 



On the other hand, the asymmetry reads 

f7r/2 fTT -, /■27r 
/O Jtt/ 



r(n)A^(n) = / '"- r IsinM^^ / / d2Q,Ar(Q,QA) 



Here 



Ma>0 



''A/i. 

A M,M' 



(30) 



E ^a:m.(KaJ'-i«-a„aJ'), 

Aa>0 



Ab = Aa - A, Al = -Aa - A and 



■n/2 r-K 



7r/2-' 



(31) 
(32) 



(33) 
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2.6 Remarks 



Two short remarks arc in order. Firstly, the expressions of the single asymmetries 
calculated in this section for particle a can be extended in a straightforward way 
to particle b. Secondly, we observe that the procedure just described, concerning 
single asymmetries, is not applicable to a decay of the type (3), but only to cases 
where both the parent resonance and at least one of the decay products are spinning 
(J, Sa(s(,) > 1/2); for example, it may be applied to 

A6^A(Ae) P{V), Ac^A P{V), (34) 

P and V denoting, respectively, a pseudoscalar and a vector particle. 



3 Double T-odd and T-even Asymmetries 



Here we define some double asymmetries, analogous to the single asymmetries of the 
previous section. As an example, we set 



A NT 





> 0) - Ar(rjv 


Tl < 0)' 




> 0) + N{T^^ 





(35) 



It is instructive to notice that the numerator of this asymmetry can be written as 

Num = [N{T^>QX>0)-NiT^>OX<0) 



\N{T^ < 0, 7f > 0) - N{T^ < 0, 7f < 0) 



(36) 



which justifies the name of "double" asymmetry. This is a T-odd asymmetry, as 
appears from eq. (35). Another asymmetry of this type can be defined as 

'n{T^ ■ > 0) - N{T^ ■ Ti^ < 0)" 



aLN 



(37) 

N{T^L ■ T,^ > 0) + N{T^^ ■ T,^ < 0)J 
Two more T-odd asymmetries can be obtained, respectively, from eqs. (35) and (37) 
by interchanging a with b. 

Viceversa, one can define five T-even asymmetries, in a quite analogous way: 
AX, Al^, A^^t, Ali' and A^^ . For example, 

'n{T^ ■ > 0) - N{T^ ■ < 0)" 



ANN _ L 
^ab — r 



N{T^ ■ Ti^ > 0) + N{T^ ■ < 0) 



(38) 
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Analytically, the first double, T-odd asymmetry reads as 



- J--K/2 A/2 



16i(-)^" 



M,M' 



The other double, T-odd asymmetry is given by 



X 



M,M' 



(39) 



(40) 



with Dh = {Xh — X'^ — l)/2. Concerning the double, T-even asymmetries, we have 



X 



X 



^JabW E « E EEQA:A,..eIV .;, ^i^;/ 



Q^AaAi^A^A' E ^Ma(^)^M'A'(^)PmM'; 



MM' 



J^jabW E « E EE01:a>„0I:a;., 

Ma, Mb Aa,A6 A^ A'^ 

:-i6) 



(41) 



n x/vx E ^MA(^^)2^^f'A'(^^)pSM'; (42) 

i'^a - A„;^A6 - A,,; M,M' 



Here 



nm'a,'{^) = 167rWj„,iy E E^<^<E^AMaM. 

Ma>0/i6>0 A 

X E 2?][;^(^^)2^^^.a(^^)pS 



'mm' 



M,M' 



Af) = E|a3[ A 51^, , 



(43) 



(44) 



11 



with Xb — Xa — A. Lastly, 



MM' 



The expressions of the asymmetries i and A^^ are obtained from, re- 
spectively, A^^ , A^^ and A^ by interchanging particle a with h. 

The asymmetries defined and elaborated in this section may be applied to a wider 
class of decays, including those of the type (3). 



4 Summary 

We have shown in the two previous sections that the expressions of the differential 
decay widths and of the T-odd and T-even asymmetries depend on interference terms 
or on the moduli squared of the "reduced" decay amplitudes. Linear combinations of 
such parameters can be obtained from experimental data, as we show in Appendix for 
several cases of interest. Indeed, we apply there the method of the moments[66, 67] 
to the distributions defined in sects. 2 and 3. As is well-known, each moment is 
factorized into a term which depends solely on the production density matrix, times 
another one which contains information only on decay amplitudes. Moreover, for each 
moment, the former factor is independent of the distribution considered; therefore the 
ratio of, say, a given moment of some asymmetry, to the corresponding moment of 
the differential decay width (24), provides linear relations among the moduli squared 
or among the interference terms of the "reduced" amplitudes. Therefore we obtain 
linear systems with respect to these quantities, and constraints to be imposed on the 
parameters which can be extracted from data. The results of the analysis that we have 
proposed can be used for determining TRV observables, as we shall show in the next 
section. In particular, the linear system obtained may be over-determined, provided 
all asymmetries defined in sects. 2 and 3 are nonzero. Unfortunately, if particle a 
or h (or both) in (4) have a strong or electromagnetic decay, some asymmetries may 
vanish because of parity conservation, as we shall see in a specific example (sect. 6). 
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However we shall establish that, in the decay considered, the nonzero distributions 
are sufficient to determining all moduli and relative phases of the amplitudes. In 
cases when this is possible, these quantities can be used as inputs for the tests of NP 
suggested in our preceding paper [49]. 



5 Determining TRV Observables 

Some TRV observables are deduced directly from experimental distributions, while 
others require some elaboration. In this section we shall consider both kinds of 
observables, starting from the former ones. 

5.1 Some TRV Observables from Distributions 

Consider the differential width of the sequential two-body decay CP-conjugated to 
(4), i. e., 

J^a b, a— )-ai 02, 6 — >■ 6i 62- (46) 
The direction corresponding to Q is 

n= (7r-^,7r + 0). (47) 

Consider the T-odd asymmetries for J. As an example, first of all, we elaborate the 
expression of the single T-odd asymmetry (fl). It is defined through 

X a\_,^ai\,^_,^ E Vi*_^imi,_^,{U)/SM'- (4J 

M,M' 



Here W is analogous to the factor (14) and p^/^/ is the spin density matrix of the 
resonance J. If the spin of the parent resonances is 0, it is natural to define a TRV 
observable as[8, 21, 25, 34, 38] 

A^^r{Q)A^{Q)+mXm, (49) 

where the -|- sign in eq. (49) is a consequence of the change of sign under CP reflection 
in the scalar product (7). For a nonzero spin of J, is a TRV observable only if 
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the production process of J is CP-conjugated to the one of J; this imphes a relation 
between the density matrices of the two resonances, practically impossible to realize. 
However, as we shall see, we can equally derive useful TRV quantities from the exper- 
imental data. To this end, we consider an ideal experiment, in which we assume that 
the production process of the parent resonance J and the decays of particles a and b 
are invariant under CP reflection. Moreover we suppose the production processes of 
J and J to be described by the same density matrix. Lastly, we take account of the 
CPT theorem, which yields W — W, and of the following properties of the V- and 
d-matrices: 

V'^.^T^i'-A^m = Vi:^{n)Vi,^,{^l), (50) 

e-\-A-.„ = -®tK,. for - K odd. (51) 

As a result, we get 

M,M' 

Substituting eqs. (26) and (52) into (49), we get 

= UfjaW E E E d!A^^MA(^)^^M'A'(^^)- (53) 
Aa,A6 A^ M,M' 

Here 

^M,M' _ (~)^° (0) rp rp _V^„SaCvSa ( rA\ 

and 

eK,\b,K = «iA6«i\ - a-Aa-Ai,«-Aj,-Ai,> (55) 
with A'„ — Aa odd. If is different from zero, one has TRV. But this imphes that 
at least one of the eA^,A6,A'„'s is non- vanishing. Although deduced under particular 
assumptions on the production process and on the decay of particle a, this condition 
depends solely on the decays of J and of its anti-particle, which are independent of the 
production process and of the successive decays of a and h. Therefore the eA„,A6,A^'s 
are TRV observables for odd — Aa. 
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One can define quite similarly the quantities A^, A^^, A^^^ and those obtained by 
interchanging, respectively, N with L and T with A^. Obviously, analogous procedures 
can be applied to these observables: A^^ gives again rise to the difference (55), while 
A™ yields 

^Xa,Xb,KK = '^Lh'^KK ~ ^-Aa-A,«-A'„-A; , (56) 
for odd A'„ — Xa and A'j, — A^. 

It is worth noting that eqs. (55) and (56) are both TRV and CP-odd, independent 
of the CPT symmetry. 

5.2 More TRV Observables 

More TRV observables can be obtained from the moduli and relative phases of the 
decay amplitudes. To this end it is convenient to pass from the helicity representation 
to the / — s one, where / is the orbital angular momentum and s the overall spin of 
the two-particle state. In particular, interesting TRV quantities may be defined as 

Here 

4 = E^lo/ CllXl ai,^ (58) 

AaAj, 

and the C's are the usual Clebsch-Gordan coefficients. The efgg,, which can be in- 
ferred from information on moduli and relative phases, generahze the TRV observable 
proposed, e. g., in refs. [8, 25]. 

5.3 Remarks 

We conclude this section with some remarks, in part connected to the CPT symmetry. 
First of all, the observables defined in this section are expressed as functions of the 

"reduced" amplitudes (16). This is not the same as obtained, e. g., in refs. [8, 24, 25], 
where the very amplitudes are involved. Therefore our analysis has picked out new 
T-odd and TRV observables, different from the analogous observables defined in those 
references. 
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Secondly, the CPT symmetry implies that the TRV observables just defined are 
of the type 



(59) 



k 



where the RkS are finite, T-even quantities and the 0jfc's are phases causing time 
reversal violation. 

Thirdly, it sometimes happens that the weak phase of a decay amplitude is neg- 
ligibly small in comparison with the strong one [48]. It is more convenient in these 
cases to consider T-odd quantities analogous to eqs. (55) or (56), but with the + 
sign in place of the - sign, or of the type (57), but with the - sign between the two 
terms. These observables are invariant under TR and, according to the CPT symme- 
try, also CP-even; they are called fake T-odd. We shall see in sect. 8 that they may 
be employed in some tests, as already proposed[38, 48]. 

Fourthly, as regards the T-even asymmetries, consider differences of the type 



where the same notations and assumptions as in eqs. (47) to (49) and (52) have been 
adopted. If significantly different than zero, these differences may consist, either of 
possible CPT-odd terms and products, or of a fake T-odd amphtude times a real TRV 
one. If CPT symmetry is assumed, these observables may be used to set constraints 
on decay models. 

Lastly, it is worth spending some words on studies of CPT violation, a very hot 
topic at present[69, 70]. Indeed, the possibihty of such a violation — connected to 
Lorentz invariance violation — was considered more than ten years ago by Coleman 
and Glashow[71, 72], who suggested experiments of neutrino oscillations. Among 
the most recent contributions, we mention the MiniBooNE[73] and Minos[74, 75] 
experiment, theoretical speculations [76, 77] and the numerous references cited therein. 
Incidentally, neutrino oscillations offer also the possibility of studying CPV and TRV 
in the leptonic sector[69, 70]. Moreover, intrinsic CPT violations, if any, have to be 
disentangled from fake violations, induced by neutrino-matter interaction [78], which 
is known as the MSW effect[79, 80]. 



A:^T{n)Al{n)-T{n)Al{n), 



(60) 
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6 An Explicit Example 



In this section we specialize the main formulae of the previous sections to the sequen- 
tial decay 

Ab (p7r~)A(//V~)j/V" (61) 

a particular case of decays (34). In this situation, the differential decay width and the 
asymmetries can be expressed as functions of the components along the unit vectors 
(6) of the polarization vector V of the resonance. Indeed, in the canonical frame, 
the Aft density matrix reads as 

= l±r.: P?i ^V.± iVy. (62) 

On the other hand, the unit vectors (6) can be expressed as functions of the angles 9 
and (f) in the canonical frame: 

Bl = {sin0cos(j),sin9sin(j),cos0), (63) 
Bt = {—sin(l),cos(l),0), (64) 
Civ = {—cosOcoscj), ~cos9sin(j), sin9). (65) 

The decay is characterized by four amplitudes, which we denote by Ai/2,1, ^-1/2,-1, 
^1/2,0 and A_i/2,oi dropping - from now on - the superscript J introduced in eq. (16). 
The differential decay width reads 

r(Q) = ^W(1 + 2VlAGl), (66) 

AGl = |«i/2,or - |a-i/2,or - |ai/2,i|^ + |a-i/2 -i|^, (67) 
Vl = P-Bl. (68) 

Now we give the expressions of some of the asymmetries: 

m)Ai{n) = ^WAa^{Bt + 2VLB^); (69) 
m)A^{Q) = iiyAa^[5R(ai/2,o«li/2,o)^iv + ^(«i/2,o«-i/2,o)^T]; (70) 
T{n)Al{n) = ^VrAa^[^(ai/2,oali/2,o)^iv-3^(«i/2,oa-i/2,o)^T]; (71) 
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T{n)Ajf{n) = -^WAa^[^{C)VN-^{C)VT]; (72) 

r{n)A^^{n) = -^wAa''[^{C)VN + ^iC)VT]. (73) 

OTT 

Here we have set 

Aa^ = lia^-a"^), (74) 

-Sl = |a;i/2,i|^ ± |ai/2,o|^ T |q;-i/2 -i^ - la-iAo^^ (75) 

Vn = V-Bn, VT = V-eT, (76) 

C = tti/2,iav2,o - «-i/2,oa-i/2 -i> (77) 



being the quantities (17) referred to the decay of the A-resonance, with positive 
(+) or negative (-) hehcity. 

The remaining asymmetries defined in the preceding sections vanish, because they 
are proportional to the quantity 

Aa^ = ^(ar-a!:i), (78) 

where V indicates the vector meson: this observable is zero, owing to parity con- 
servation in strong and electromagnetic decays. However, information that we can 
extract from the non-vanishing asymmetries is sufficient to determining the moduli 
and relative phases of the reduced amphtudes. Indeed, the moduh can be inferred 
from T{Q) and Aj(^{Q); moreover the relative phase between Q!i/2,o and '^^^ ^le 

extracted from (^) or from A'j^{fl); lastly, the sine and cosine of the other two 
relative phases are related, respectively, to Q{C) and to 3?(C), whence the two phases 
can be deduced. 

7 Relations to Previous Works 

In this section we compare our results with those of other authors. In particular, 
as we shall see, our choice of considering correlations constructed exclusively from 
momenta results to be rather general. Indeed, such correlations turn out to be spin 
dependent, yielding results equivalent to those which involve one or more spins. A 
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strong indication in this sense comes from the observation that the asymmetries 
defined in sects. 2 and 3 vanish if both the parent resonance and the decay products 
are spinless. 



7.1 Hadronic Decays 

Let us focus on decays of the type (4), such that J and or/and Sf, are nonzero. For 
example, consider decays (34), studied by several authors[24, 40, 50, 51, 55] and also 
detected experiment ally [81, 82, 83, 84]. We define the following products: 

Sa(6) • Gat, Sa(b) ' ^T: Sa(b) ' ^L, (79) 

whose mean values are the components of the polarization vector of one of the two 
decay products, a or 6; we denote such components, respectively, as Pjv, Pt and Pl- 
Pjv is T-odd, while Pt and P^ are T-even. If both a and h are spinning, one can also 
define polarization correlations, like, e. g., 

Ptn ^ {Sa-er Sb-BN)- (80) 

This is the case of the decays (34), if a vector meson V is involved. Here the normal 
component Pn of the polarization of the final fermion reads as [61] 

m)Pi:^{n) = i-(G^ + AG%Vl), (81) 

where 

G% = 3fj(Ai/2,iA:^i/2,o + ^-1/2,-1^1/2,0)' (82) 

AG^ = -2^{A^/2,^AU/2,o- A-i/2,-iAl/2,o). (83) 

As regards the polarization correlation Ptn, we have[61] 

r{n)PTNm = -^{AGtn + GtnVl), (84) 



with 



Gtn = 2$5(A_i/2,_i/lt/2,o + ^1/2.1^-1/2,0), (85) 
AGtn = ^^(A_i/2,-iAt/2_o-A/2,i^li/2,o)- (86) 
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As can be seen, such polarization formulae consist of linear combinations of the pa- 
rameters which appear in the single and double asymmetries defined in sects. 2 and 3. 
Also the other components of the polarization vector and the polarization correlations 
are related to our present results in a similar way[61]: see also sect. 6. It is reasonable 
to think that this result is not peculiar to the decay considered; we conjecture that 
the asymmetries defined in the present paper provide, in principle, information on 
the components of the polarization vector and polarization correlations. Even, the 
results of our analysis may be quite suitable for determining the components of the 
polarization vector of a decay product which carries a spin greater than 1/2; indeed, 
it is not so easy to measure directly this polarization. 

7.2 V1V2 Decays 

As we have seen in sect. 3, double asymmetries include decays of the type (3), for 
which our analysis yields the following T-odd observables: 

^(cKiiCKoo) ^(q!_i_iq;qq). (87) 

A linear combination of these terms corresponds to the usual T-odd quantity[21, 25, 
34, 38, 48] 

^{A^AD, (88) 

where the A's are the decay amplitudes in the transversity representation. In fact, 
we have the following relations: 

Aq = Faoo, ^11 = --^F(q;ii q;_i_i) A±_ = - «-i-i), (89) 

with 

F2 = ^ \Ax\^ = |Ao|2 + + \A^\\ (90) 

A 

Moreover, for the decays considered, eq. (56) yields the following TRV observables: 

^(aiiaJo - ^-1-1^00) and S(a_i_iaJo -^11^00)5 (91) 

which are linearly related to '^s^A^Aq), to Q'(74||74q) and to their CP-conjugated quan- 
tities. 
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Some remarks are in order. 

a) It is worth stressing that also in this case our method yields results which turn 
out to be very similar to those found by the authors just mentioned, although they 
adopt the T-odd product 

Pa • s„ X Sb. (92) 

b) The usual observables /l and /t[21, 23, 35, 38, 39, 48], concerning respectively 
longitudinal and transverse polarization of the vector mesons in a decay (3), and used 
for testing the SM[38, 39, 47], are simply related to our "reduced" amplitudes: 

Jl ^ \oioof, h = |q;ii|^ + \a-i-i\^. (93) 

c) The T-odd observable[25] $>(Aj_Ap does not appear among the parameters 
of our analysis. This is due to the fact that we have integrated over some angular 
variables in the sequential decay, whereas one usually adopts a full angular analysis [16, 
17, 28, 29, 85], in terms of two polar angles and an azimuthal one [42, 44, 86]. But if, 
as discussed in sects. 4 and 6, our method allows to extract the moduli and relative 
phases of all of the decay amplitudes, the above mentioned term can be deduced 
indirectly, as well as those which do not appear as parameters of our analysis. Our 
procedure may be suitably followed in cases where the statistics is not so rich as in 
Babar[28] and Belle[29] experiments. 

d) The TRV observables proposed by previous authors, e. g., 

'^{A^Al + A^Al), (94) 

involve just the decay amplitudes. But we observe that the proportionality constant 
F relative to the CP-conjugated decay is different than F . This implies that the 
observable (94) is substantially different than (91). More generally, it confirms once 
more that the TRV observables defined in subsects. 5.1 and 5.2 are really different 
than those proposed in the literature. 

7.3 Detection of Top Decays and New Particle Decays 

Also some of the top decays[9, 58, 59] and of the possible new particle sequential 
decays[18, 20, 22, 60] are studied by means of triple products. In particular, it is 

21 



pointed out that correlations from three or more momenta[18] appear to be more 
appropriate for experimental detection and that the maximal asymmetry is obtained 
by operating in the rest frame of the decaying particle[18]. 

8 Tests of Self-Consistency and of SM 

Here we propose some tests, based on the observables illustrated in sections 2, 3 and 5. 
Moreover, as shown in sect. 4 and in Appendix, the analysis that we have illustrated 
allows to determine quantities of the type 

\A^\\ ^{A^AD, ^{A,Al), (95) 

where Ai and A2 are two different helicity amplitudes of a given decay mode. Now 
we propose various kinds of tests, based on the knowledge of such parameters and of 
others, related to them. 

8.1 Self-consistency Tests 

A first kind of tests is a self-consistency one. We suggest to exploit the identities of 
the type 

[^{A^ADY + ["^{A.ADf = \A,\^\A2\\ (96) 
as constraints on the parameters to be extracted from data. 

8.2 Tests of Moduli and Phases 

As told in the introduction, it is always convenient to look for observables whose 
values are small according to the SM, since they may give a clear indication of NP. In 
particular, it appears suitable to study decays such that the SM predicts neghgibly 
small TRV or/and direct CPV. 

Some tests in this sense were akeady proposed in a preceding paper [49]. For 
instance, we suggested to consider the asymmetries 

^ |a^a^^^-a^ and = 14^^^4^-^=^^^- (97) 

$A„A, + ^-A„-A, I^AaAj' + |^-A„-aJ' 
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Here $A„Ai, is the relative phase of the amphtude ^a„A6 to a fixed one, taken as a 
reference. 

Aside from that, as aheady observed by various authors (see, e. g., ref. [8, 21, 24, 
25, 34]), a hierarchy among the CP-odd or T-odd observables has to be cstabhshed, 
according to their sensitivity to NP. For example, consider the interference term 
A1A2. In order to evaluate qualitatively the behavior of this term, we assume a quite 
simplified model, such that each amplitude is of the form 

A = Te*^^ + Pe'^'^p+^^l (98) 

Here we have omitted, for the sake of simplicity, the index 1 or 2; moreover the former 
term corresponds to the tree contribution and the latter to the penguin graph; in 
particular, the cu's (assumed to be helicity independent) and the 5's are, respectively, 
the weak and the strong phases. As usual, we assume Sp to he attributed exclusively 
the absorptive part of the penguin diagram[7]. Now we substitute eq. (98) into the 
interference term and compare it with its CP-conjugated quantity. As a result we get 

^{AiAl-AiT^} = -2sinAuj{\Pi\\T2\sinSpi + \P2\\Ti\sinSp2), (99) 
•^{AiAl-AiT^) = 2sinAuj{\Pi\\T2\cosSpi-\P2\\Ti\cosSp2), (100) 

where Au — OJp — OUx. Then we conclude that, if the strong phase is small, the real 
part of the interference term is surely negligible, while the imaginary part may be 
sizeable, provided the difference ||Pi||T2| — |P2||rL|| is sufficiently large[7]. Therefore 
the imaginary parts of the complex quantities (55) and (56), as well as the observable 
(57), appear favourite in the search for possible clues of NP. 

8.3 Helicity tests 

Here we suggest an interesting test of the SM, as a generalization of the one proposed 
in ref. [48]. This test is not connected to interference terms, but it is an important 
byproduct of the helicity representation, used in our treatment. Indeed, the SM pre- 
dicts positive helicity amplitudes to be strongly suppressed with respect to negative 
helicity ones, if heavy quarks are involved in a decay. This is verified at an empirical 
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level, although theoretical arguments are not so sound[39]. Then violation of this in- 
equality may be an indication of NP. Incidentally, if the decay products are spinning, 
and at least one of them is a vector meson, it is interesting to compare the 0-helicity 
amplitudes with the corresponding ones with ±l-helicity amplitudes. Factorization, 
generally satisfied by the tree contribution, predicts that the 0-helicity amplitudes 
are much greater[39]. On the contrary, the penguin contribution, for which factoriza- 
tion fails, yields comparable 0- and ±l-helicity amphtudes[39]. Therefore the ratios 
between such amplitudes indicate the relative weight of penguin to tree term. It is 
worth noting that such kinds of tests can be suitably performed by means of the fake 
T-odd observables introduced in sect. 5, especially if the weak phases are neghgibly 
small[48]. 

8.4 and Ac Decays 

It is useful to consider the application of the previous tests to decays of heavy baryons, 
for example to those described by eq. (34), in view of the forthcoming LHCb data. 
Contributions in this sense, especially if one of the decay products is a vector meson 
V, have been already given[24, 40, 49, 50, 51, 55, 61, 87]. According to the helicity 
tests, we expect |Ai/2,i| << |^-i/2,-i| and |^i/2,o| << |^-i/2,o|) the A^s being the 
decay amplitudes introduced in sect. 6. 

Moreover the ratios |A_i/2,-i/^-i/2,o| and |^i/2,i/^i/2,o| may give indications on 
the proportions with which penguin and tree diagram contribute to the decay ampli- 
tudes. Let us consider the case oi V — J/'0[61]. Assuming amplitudes of the type 
(98), the SM predicts 

Te''^^ = HcKT, Pe*(-p+5p) ^ ^ (101) 

Q 

Here the V^s are elements of the CKM matrix, T, P and T' are real positive numbers 
and q = u,c,t. Taking into account the orthogonality condition 

and the fact that VbuV*^ is negligible in comparison with the other two terms, we 
conclude that the phase difference Acu is quite small according to the SM. However, 
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recent experimental results, concerning both Bg — Bs mixing[30, 31] and direct CP 
violation asymmetry in S — > Xtt decays[27, 28, 29], indicate that this phase could 
be considerably larger. Therefore it is worth applying our analyses to these decay 
modes. Incidentally, we signal that a commonly used model[88, 89, 51] predicts P/T 
- 0.134[51]. 

Even more intriguing is the case oiV — p°, a;[87], or the decay Att". Indeed, 

here the penguin contribution is greater than the tree; therefore the interference 
between the two terms, which includes the CP-violating phase, is quite relevant. 
Moreover, the diagrams involved are quite similar to those which occur in the decays 
B — > ir7r[27, 28, 29]. Therefore it is worth doing efforts for revealing such decay 
modes of A?, and for applying to them the tests suggested in this section. 

As a further example of decay, consider [90] 

Ac A7r+. (103) 

In this case no CP violation has been observed [90]; however an analysis of the triple 
product asymmetry is suggested, since one has to do with a decay where the relative 
strong phase of the two amplitudes is quite small. Obviously, our method appears to 
be appropriate also in this case. 

9 Conclusions 

We have elaborated some methods for analyzing hadronic sequential two-body decays 
involving more spinning particles. In particular, we have suggested to investigate sev- 
eral distributions, based on T-odd and T-even, simple or double correlations. Unhke 
other authors[43], who try to disentangle different CP eigenstates, we exploit just 
the interference between such eigenstates. The decays considered offer a richer range 
of observables sensitive to CPV and TRV. They may also help to find hints to NP, 
provided we focus preferably on observables for which the SM predicts quite small 
values. Now we exhibit the highlights of our analysis. 

a) Our main result is that, given a set of data concerning the above mentioned 
decays, one can always infer a set of TRV observables, even after integrating over 
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all variables, except for the direction of one of the two decay products in the initial 
decay. We stress that this would not be possible if less than two of the particles 
involved in this decay arc spinning. Some of these observables are especially sensitive 
to TRV and to possible clues beyond the SM. Among them we signal the imaginary 
parts of the interference terms like (55) and (56), which are rather small according 
to the SM predictions. In particular, the study of the interference terms in the decay 
modes A;, — >■ Ap^,uj appears very interesting in this sense. 

b) We have shown that the T-odd correlations based just on momenta (three or 
more), which may be more easily adopted in an experiment, provide the same results 
as those involving one or more spins. 

c) Thirdly, T-even observables - quite helpful, although a bit neglected in the 
current literature - may be combined with their CP-conjugated ones, so as to in- 
clude either possible CPT-odd terms or products of a fake T-odd amplitude, times 
a real TRV one. These observables, as well as fake T-odd ones, may be used to set 
constraints on the T-odd terms caused by strong interactions, typically spin-orbit 
ones. 

d) Our treatment recovers, as particular cases, some methods suggested by other 
authors. For example, the analysis of the helicity amplitudes, rather efficient as a 
test of the SM, is an important byproduct of our choice of the helicity representation. 
Furthermore, new TRV observables have been picked out, the "reduced" amplitudes, 
which generalize the polarization parameters and fx, used for vector mesons. 

e) Lastly, our method may allow, at least in some cases, to determine the moduli 
and the relative phases of the decay amplitudes. These quantities allow, in turn, to 
infer further, especially sensitive TRV observables and the polarization vectors of the 
decay products, otherwise difficult to determine for, say, a vector meson. 

The analysis proposed here could be usefully apphed to forthcoming data, espe- 
cially at LHCb, but also in experiments where the available statistics is not so rich 
and abundant. 
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Appendix 

We describe a method for extracting from data the moduh squared of the hehcity 
decay ampUtudes and some interference terms between them. The method is based 
on the moment expansion of the various distributions defined in the text, that is, 
r(17), (f2)r(f2), etc., which we denote generically by F[Q). Each such distribution 
contains the product T^ma{^)T^m'a'{^) of Wigner "D-matrices. Therefore our starting 
point is the well-known relation 

'^Ma(^)^M'A'(^) = Cjf/ j^M VA ^N-t/ (^) ■ (A. 1) 

L 

Here the C's are the usual Clebsch-Gordan coefficients. Inserting this into F{Q), we 
get an expansion of the type 

F{Q) = Y: Hi^MUQ). (A. 2) 

Here 

ttJ — fJ* fJ fJ* _ „(0) (~iJ L J ( k 

MM' 

and the //,^'s vary from distribution to distribution, as we shall specify below. They 
depend only on the decay amplitudes, while the coefficients depend only on the 
production reaction of the parent resonance. 

A.l - Extracting Moduli of Decay Amplitudes 

Now we give the expressions of the coefficients //^ for three distributions, which 
depend solely on the moduli of the decay amplitudes, that is, r(r2), T {VL) A^^p'^{VL) and 
V{yL)A^^{yL), whose expressions are given, respectively, by eqs. (24), (30) and (43). 
We have 

fi^P = MjW5.o E \<x,?Ci^L (A. 4) 
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Here the various symbols introduced are defined in the text. We just reproduce those 

A6- 



expressions which contain the "reduced" amphtudes a'l 



= ^ E ^t,S\<x,? - l«-A„A;r), (A. 7) 

^ A„>0 

^A/!aMb ^ E I^AaAb^^AaMa^AlMb- (A- 8) 

A„ 

Moreover the upper indices of the coefficients //^ refer to the specific distributions 
considered; in particular, the index £ refers to "longitudinal asymmetry", not to be 
confused with the order L of the moment. Lastly, the coefficient /^^ can be defined 
analogously to eq. (A. 5), by interchanging index a with index h. Now we consider 
the ratios 

rf = Hi^,f/HEl (A. 9) 

r'l = Hi^^!^/Hi^Z (A. 10) 

= Hi^^P/Hi^P,. (A. 11) 

Here, again, the upper indices of the H'^^s arc typical of the distribution considered. 
First of all, it is important to note that the ratios (A. 9) to (A. 11) arc independent 
of N . This constitutes a check for the moments of the distributions. Furthermore, 
by recalling the relations above, we obtain the following linear system in the moduli 
squared of the various ocj^^)^^'- 

[.xfci^i - ^^-^T.^<i:^<ci^L (A. 12) 

:xfci^l = ^^^T.^<Y.^<fi^L (A. 13) 



' L 


E 




Aa,A(, 




E 






U 


E 




Aa,Aj, 



•'^•^ (Ua>0Al6>0 A 

The parameters Aa^^^'j,,, relative to the secondary decays, are generally known. The 
linear system (A. 12)-(A. 14) is not homogeneous, owing to the constraint 

E KaJ^ = 1- (A. 15) 
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It is generally over-determined, as we shall see in some examples of interest, provided 
all of the asymmetries are nontrivial. 

A. 2 - Inferring Relative Phases of Decay Amplitudes 

The relative phases of the decay amplitudes may be determined starting from 
the other distributions considered in the text, where interference terms between am- 
plitudes with different helicities are involved. Here we limit ourselves to two of the 
T-odd distributions, but the method we suggest can be extended to any distribution 
of sects. 2 and 3. 

According to our formalism, the moments related to the distribution r{0,)A-^{fl) 
(see the first eq. (A. 3)) yield the following decay coefficients //j,: 



X 



\a Xb K K ^ '^JK'^b '^bl 



X <xy>lK^ (A- 16) 

A = K-\, A' = A;-A;, L>„ = (A„ - a; - 1)/2. (A. 17) 

Here, as in the text, the symbol means that Aq — A^ or A^ — A|, is an odd quantity. 
Analogously, the distribution T(p)A^^ {Q) yields 



IJ-b 



X 



Xa Xb A^ Afe-Aft 



X <Xb<x',- (A. 18) 

Moreover — (A^ — A'j, — 1)/2 and eqs. (A. 17) hold. Other equations can be obtained 
from the other distributions considered in sects. 2 and 3. Once the moduli of the 
decay amphtudes are known, the ratios Hj}^J^ / Hj}pQ, Hjj^^^/Hj}pQ and the other 
ratios obtained from the various distributions depend only on the relative phases 
of such amplitudes. It can be checked that, in the cases of interest (see the next 
subsection), a complete set of n — 1 independent relative phases - n being the number 
of helicity amplitudes - is present in the equations that we obtain from these ratios. 
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Again, if the asymmetries are all nontrivial, this system is over-determined in the 
cases of interest. 

A. 3 - Some Examples 

Wc present some cases of interest of two-body decays in Table 1, where we give the 
number of unknowns and of available equations. As regards relative phases, we have 
taken into account only the conditions deriving from the knowledge of the moments 
defined in subsect. A. 2. 

Table 1: Moduli and relative phases in different channels 

Decay N^n N; 

1^0 110 

0^11 3 4 2 2 

0^ 1/2 1/2 2 4 1 2 

1/2^ 1/2 2 7 1 4 

1/2^ 1/2 1 4 7 3 4 

3/2^1/2 1 6 13 5 8 

Here N^(j,) is the number of moduli (relative phases) of of the decay amplitudes 
in the various channels and N^^^^ the respective numbers of equations. 
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